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We use the Interacting Instanton Liquid Model (IILM) as a tool to study the role played by the 
chiral interactions in the lowest-lying vector and axial vector meson resonances. We find that narrow 
a\ and p meson resonances can be generated by instanton-induced chiral forces, even in the absence 
of confinement. In the IILM, these hadrons are found to have masses only about 30% larger than the 
experimental value and small width < 10 — 50 MeV. This result suggests that chiral interactions are 
very important in these systems and provide most of their mass. We explore the decaying patterns 
of the p meson, in the absence of confinement. We argue that, in our model where only chiral forces 
are switched on, this meson decays dissociating into its quark anti-quark constituents. 



I. INTRODUCTION 



o 
o 

a 

)—} , Hadron spectroscopy provides basic constrains on the structure of non-perturbative QCD dynamics. Since gluons 

■ interact with quarks in a way which depends dramatically on the quark mass, we expect the light- and heavy-hadron 
£SJ , sectors of the spectrum to be sensitive to different dynamical correlations. For example, while charmonium spectrum 

is well understood in terms of linearly rising confining potential and perturbative gluon exchange, the description of 
the light-hadron spectra in terms of QCD degrees of freedom must take into account also the interactions responsible 
for spontaneous breaking of chiral symmetry. In particular, we expect the mass and the structure of lowest-lying 
hadrons such as pion, nucleon, vector and axial-vector mesons to be strongly influenced by the chiral interactions, 
. because the splitting between parity partners in this part of the spectrum is as large as 500 — 600 MeV. On the other 

■ hand, Regge trajectories suggest that resonances with large angular momentum are predominantly influenced by the 
physics of color confinement, and it has been conjectured that chiral symmetry breaking may be even restored up in 

t — ■ the spectrum (see e.g. [1| and references therein). In this paper we address the following questions: are any of the 
light hadron resonances generated predominantly by chiral forces, with confinement playing a sub-leading role? Could 

■ any of such hadrons exist even in the complete absence of confinement? In this case, would such systems still decay 
rS ! predominantly in colorless hadrons? What is the microscopic dynamical mechanism underlying the splitting between 

' , the different lowest-lying multiplets? 

The answers to these questions root in the non-pertubative sector of QCD. Although lattice field theory provides 
the only available ab-initio tool for computing non-perturbative QCD correlation functions, the mechanism by which 
the hadron structure arises is not directly evident in such a framework. In particular, it is difficult to disentangle the 
contribution to the correlation functions arising from the different types of interaction, such as the confining forces, 
the chiral forces, and the perturbative gluon exchange. Hence, in order to gain physics insight, in the present work 
we focus on the role played by chiral symmetry breaking and we do so by restricting the path integral to a sum over 
gauge field configurations which generate the near zero-mode zone of the Dirac operator. To this end, we adopt the 
interacting instanton liquid model, developed in Instantons have been long argued to be the dominant fluctuations 
generating the zero- mode zone of the Dirac operator, hence providing the correlations which break spontaneously chiral 
symmetry. In the IILM, the QCD path integral over all gluon configurations is replaced by an effective theory in 
which the gauge fields accounted for are those generated by integrating over the positions, color orientations, and 
sizes of singular-gauge instantons. 

Several features of the instanton picture have been observed in a number of lattice studies 

HHHH- In particular, 

it has been observed that chiral symmetry breaking in QCD is strongly correlated with smooth lumps of topological 
gauge fields, with a shape compatible with that of singular-gauge instantons 

Instanton-induced correlations in hadrons have also been studied through several phenomenological model calcula- 
tions. It has been shown that instanton models provide a good description of the mass and the electro- weak structure 
of pions, nucleons and hyperons [H, EH EH EH EH EH El, EE EH EH ■ In a recent work [l9| , we have used the IILM 
to study the dependence of the Dirac spectrum and of the pion and nucleon masses on the current quark mass. We 
have found that the dependence of these observables on the pion mass is in quantitative agreement with both chiral 
perturbation theory and lattice simulations. 

The main shortcoming of the instanton liquid model is that it does not lead to the correct behavior of the Wilson 
loop, i.e. it does not predict a linearly rising potential between static color sources, at large distances. This fact 
implies that, while singular-gauge instantons accounts well for the non-perturbative QCD correlations at the chiral 
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symmetry breaking scale A x ~ 1 GeV, they do not generate sufficient long-range correlations at the confining scale 
~ ^qcd ~ 0.2 GeV. As a result, we expect the chiral forces generated by these topological fluctuations to play 
only marginal role in heavy-quark systems — which are insensitive to chiral dynamics — or in highly excited hadrons 
— whose wave functions are expected to extend for several fm — . On the other hand, just because of the lack of 
confinement, the IILM can be used as a tool to single-out the contribution of chiral forces in the different hadronic 
systems. By studying where the model works well, we can identify the properties of the hadrons which are almost 
completely determined by the chiral dynamics. Conversely, by studying where the model fails, one can in principle gain 
information about the structure of the additional interaction which is needed in order to reproduce the experiment. 
In this context, a\ and p resonances represent the ideal test systems: being in between the pion — which is strongly 
related to chiral symmetry breaking — and highly excited states — which are dominated by color confinement — they 
are expected to be sensitive to both chiral and confining interactions. In addition the splitting between these two 
resonance masses is a direct measure of the effect of spontaneous chiral symmetry breaking. 

An exploratory study of the contribution of instanton forces to the lowest- lying hadrons was performed in [l(| [H| , 
where several hadronic Euclidean point-to-point correlation functions where calculated in the ILM up to sizes of the 
order of several fractions of a fm. It was shown that instanton-induced chiral forces are very strongly attractive in the 
nucleon and pion, but much weaker in the p and A. Their study provides indications that hadron resonances might 
exist in the instanton vacuum. In fact p, a\ and A masses extracted from a fit of the corresponding point-to-point 
correlation functions using a pole-plus-continuum ansatz for the spectral function, were found to be within 20 — 30 % 
from their experimental value. 

In this model, the dynamical origin of the splitting between the hadron multiplcts which are not connected via chiral 
transformations is well understood in terms of the diluteness of the instanton liquid. In fact, the instanton zero-mode 
contribution to the two-point correlation functions associated to the pion and to the nucleon comes at the lowest-order 
in the instanton liquid diluteness k ~ 10 _1 . On the other hand, the analog contribution to the correlation functions 
associated to vector- and axial-vector mesons and to decuplet baryons masses comes at the next-to-leading order, i.e. 
to 0(k 2 ). The splitting between chiral partners is provided by the spontaneous chiral symmetry associated with the 
de-localization of the fermion zero- modes (see [2(| and references therein). 

In this context an interesting question is whether, in the presence of only instanton-induced chiral forces, the a\ and 
p resonances predominantly tend to decay into colorless pions, or rather dissociate into free quarks and anti-quarks. 
In the present work we address this and other related issues concerning the contribution of chiral dynamics to the 
structure of hadronic resonances by computing momentum projected correlation functions at Euclidean times up to 
~ 1.2 fm and we focus on the effective mass plot. 

This method is usually employed in lattice QCD simulations to extract the masses of the stable bound-states only, 
i.e. when the lowest-lying hadron contributes to the spectral density through a (^-function. 

In this case, the effective mass plot displays a flat plateau in the large Euclidean time limit. We extend this 
method to include the case in which the lowest-lying hadron is not a stable bound-state but rather a resonance with a 
finite width. In this case then effective mass plot displays an approximatively linear and mild fall-off at intermediate 
Euclidean times. This behavior can be distinguished both from the flat plateau associated to stable bound-states and 
from the strong exponential decay associated to the filtering of a perturbative continuum. As a result both the mass 
and the width of the lowest-lying mesons and baryons can be extracted. 

In addition, some information about the decay patterns of these resonances can be gained by studying the behavior 
of the effective mass in the asymptotically large Euclidean time. In this regime, the effective mass converges to the 
smallest eigenvalue in the transfer matrix. If the simulation box is sufficiently large, the lowest eigenvalue of the 
transfer matrix in the p and a\ channels corresponds to the invariant mass of the branch-cut singularity in their two- 
point function. In QCD these are expected to be at the threshold for decaying into two- and three- pions, respectively. 
In models which do not account for confinement, these hadrons could also decay into free quarks and the branch-cut 
may be shifted. 

We apply the effective mass plot analysis to determine the mass of the lowest-lying vector and axial-vector meson, 
using several values of the quark mass. Our results clearly show that the instanton-induced chiral forces are sufficiently 
strong to generate the p meson and the a\ meson resonances even in absence of confinement. 

On the other hand, while our previous studies have shown that the nucleon and pion masses are correctly reproduced 
in this model [19 ]. the vector and axial vector mesons resonances are found to be about 30% heavier than their 
experimental value, with M p ~ 1 GeV and M ai ~ 1.7 GeV. These results support a picture in which the mass and 
binding of the nucleon and pion is almost entirely provided by chiral forces, while the mass of the vector and axial 
mesons receive a significant contribution from confinement interactions of the order of 30%. 

For some of the quark masses we have used, the resulting mass of the p and a\ resonances turn out to be significantly 
smaller than the threshold for decaying into three- and two- pions, respectively. In the presence of confinement such 
hadrons would be stable. However, in our IILM calculations we find a p meson width of 10 MeV. Our interpretation 
of these results is that, in the presence of chiral interaction only, the p meson eventually tends to dissociate into its 
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FIG. 1: Left panel: parametrization of the ALEPH Collaboration [2l[ data for the vector spectral density. Right panel: the 
corresponding effective mass plot 

constituents. 

In the next session we present our phenomenological approach to extract information about hadron resonances from 
effective mass plot. In session HTT1 we present and discuss the results of our calculations. Session HVl is devoted to the 
summary and conclusions. 



II. RESONANCES AND THE EFFECTIVE MASS PLOT 



In QCD the information about the hadron spectrum is encoded in the two-point correlation functions, defined as 

n fl (x,T) = (o|t[j h (x,t)4(o,o)]|o>. (i) 

</# (x, t) is an overlap operator that creates states with the quantum numbers of the hadron H . The lowest-dimensional 
overlap operators generating states with quantum numbers of 7r, p, a\ mesons are: 



Jh{x) = q{x) T H q{x), 
In the following we consider the effective mass M e J f (r), defined as 



M e H ff (r) 



lim — — In 

At->o At 



G h (t) 



_G h (t + At) 

where Gh (t) is the zero- momentum-projected hadronic two-point function, 

G h {t) = J d 3 x n H (x,r), 
which can be written in the spectral representation: 

d.s 



Gh(t) 



p H (s)e 



(2) 
(3) 



(4) 



(5) 



(6) 



where Ph(s) is the spectral function. 

In the large Euclidean time limit, the effective mass filters-out the lowest singularity in the two-point function, i.e. 
the smallest eigenvalue of the transfer matrix. If the lowest-lying state in a given channel is a stable hadron, then the 
two-point function develops a pole at the bound state mass below the threshold of the branch-cut associated to multi- 
particle production. Then, the effective mass asymptotically approaches the value of the mass of the bound-state: 



lim M e J f (r) = M H 



(7) 



In QCD the only bound states are pions and nucleons. In all the other channels, the two-pint functions displays 
branch-cut singularities only. 
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It is instructive to study the behavior of the effective mass plot if the spectral function ph(s) displays a narrow 
resonance with a finite-width, emerging above a continuum back-ground at small s and converging to the asymptotic 
perturbative continuum, in the large s limit. As a working example, we consider the effective mass plot for the vector 
meson channel. In this case, the spectral function can be extracted from the ALEPH Collaboration data [2lJ for t 
decays in two pions. A reasonable parametrization of such data can be constructed from a Breit-Wigner function for 
the p-meson resonance, supplemented by a term simulating the perturbative continuum (see Fig[TJ left panel) [22I ]: 

n {<,) = r p ( r p/ 2 ) 2 1 9i (a\ 

Pp[ '> 0l (r p /2)2 + (Vi-m p )2 + I + exp[(£ - V^)/0.2] W 

The right panel shows the effective mass plot obtained from a phenomenological parametrization of the two-point 
function, using Eq.s (jlj, ([6]) and (jHJ. At small Euclidean times, r < 0.4 fm, the effective mass M^Jt) drops 
exponentially. Such a rapid fall-off is due to the exponential suppression of the perturbative continuum of excitations 
induced by the propagation in the imaginary time. At larger Euclidean times, 0.5 fm < t < 2 fm, the effective 
mass displays a linear, nearly flat region. In this regime, the spectral representation of M^j(t) is dominated by the 
p-meson resonance peak. In fact, it is easy to check that, in the limit of vanishing width, one recovers a completely 
flat straight line, i.e. the familiar signature of a stable bound-state. Eventually, at even larger r, the effective mass 
slowly converges to the threshold energy for multi-particle production 1 . We stress the fact that the effective mass 
plot analysis is much more efficient than the corresponding point-to-point correlation function study in distinguishing 
a resonance peak from a stable bound-state. 

From this discussion it follows that it is in principle possible to extract the width of the resonance from the effective 
mass plot in the intermediate Euclidean time region. We note that, in (unqucnchcd) lattice simulations, the stability 
of all hadrons except nucleon and pions depends on the size of the simulation box and on the value of the pion mass 
used. For example, at large pion masses > m p the p-meson is a stable state, because there is no phase-space 
available for decaying into two pions. In this case, one can simply read-off its mass from the plateau in the effective 
mass plot, at large Euclidean times. 

On the other hand, for sufficiently small pion masses, the phase-space for decaying into two-pions opens up and the 
p-meson appears in the spectral density as a resonance. In this case, the smallest eigenvalue of the transfer matrix 
filtered-out by the propagation in imaginary time is related to the two-pion p-wave state. Note however that, if the 
periodic box is too small, the quantization of momentum may shift the p-wave two-pion state threshold above the 
p-meson mass. As an example, let's consider a simulation performed in a box with size of L = 2.5 fm and with a pion 
mass of 500 MeV. In this case, the smallest non- vanishing unit of momentum is 2ir/L ~ 500 MeV, and the threshold 
for decay into a two-pion p-wave state is at + (ir/L) 2 ~ 1.1 GeV which can be above the p-meson mass. 

Note that in our IILM calculations we do not have to worry about effects related to quantization of momentum, as 
we do not adopt periodic boundary conditions. Instead, we choose simulation boxes which are sufficiently large for 
the integrand Hh (x, t) in the momentum projection integral |5j) to become very small and negligible near the borders 
of the box 2 . Under such conditions, the lowest point in the branch-cut for the p meson two-point correlation function 
is located at the threshold for two-pion production, i.e. -,/s = 2M W . 

Let us now discuss the axial-vector channel. In this case, the hadronic current has an overlap with both the 
pion state and the a\ resonance. A rough parametrization of the ALEPH Collaboration data [2l[ for r decays into 
three-pions (see the left panel of Fig. [5]) leads to the spectral function [22| : 

P {s) = c? jJtfnr V2 ~ fi m W a ~ m " ) + tt — u? rwn9v ( 9 ) 

(T ai /2y + (V s - m ai y 1 + exp^o - VSJ/0.2J 

where the pion pole arises from the matrix element (0\j£(0)\n) = ip^f-n- We note that the pion contribution to this 
spectral function comes with an opposite sign with respect to that of the a\ resonance. Using ^ and dialing the 
physical value for and /„•, we obtain the effective mass plot shown in the central panel of Fig[2j This plot displays 
a structure which is qualitatively similar to that of the vector-meson channel. On the other hand, dramatic differences 
emerge when the width of the a\ is reduced (for instance because of heavy pions in the spectrum). For example, if 
we reduce r oi by one order of magnitude we obtain the effective mass plot shown in the right panel of FigO In the 



Note that, since our simple phenomenological parametrization JSJ does not vanish below the two-pion threshold, s = \ / 2m-^, the 
resulting effective mass converges to in the asymptotically large Euclidean times. 
2 Note that this is different from imposing Dirichlet boundary conditions, as we do not impose wave-functions or correlators to vanish at 
the border of the box. 



FIG. 2: Left panel: parametrization of the ALEPH Collaboration data [5l|] for r decays in three pions using Eq. (J9j . Center 
panel: the corresponding phenomenological effective-mass plot for the ai meson. Right panel: effective-mass plot for the a± 
meson with reduced width F ai = 0.03 GeV 



limit of stable ai (vanishing width), a singularity develops at the Euclidean time 



M ai - NU 



log 



A' 

ai 



(10) 



where A ai is the coupling of the axial-vector current to the a\ state. This is a consequence of the cancellation between 
the contributions of the pion and axial- vector poles in the denominator of Eq. (|4|) . 

In the next session we compare these phenomenological representations of the effective mass plot with the results 
of calculations performed in the IILM. 



III. RESULTS AND DISCUSSION 



In the IILM, hadronic correlation functions are evaluated by means of Monte Carlo averages over instanton ensamble 
configurations. The only phenomenological parameters of the model are the instanton average size p = 0.33 fm and 
the dimensionless strength of the instanton-antiinstanton bosonic short-distance repulsion (for a concise review of 
this model see [l9|], for an extended treatment see [13] )■ In the present calculations we used five sets of ensemble 
configurations, corresponding to quark masses ranging from 20 to 90 MeV and we estimated statistical errors using 
jackknife technique, with bin size of 10 configurations. In order to isolate the instanton-induced chiral interactions, 
we have adopted the so-called zero-mode approximation, in which the part of the quark propagator which does not 
receive contribution from the instanton zero-modes has been replaced by a free propagator (for further details, see 
the discussions in [20 ] and [19| ). For comparison, we have performed the same calculations also including the non- 
zero-mode part of the propagator and we have not found significant differences, hence instanton-induced correlations 
not associated to the chiral zero-mode zone play only a marginal role. 

Our results for the calculation of the p meson effective mass for several quark masses are shown in Fig[3J where 
they are compared with the best fit obtained from the phenomenological representation of the spectral function (jHJ). 
The IILM points can be very well interpolated throughout the entire region 0.4 fm < r < 1.2 fm. On the other hand, 
it should be noted that a fit of these points using a spectral function which does not account for a narrow meson 
resonance would be inconsistent with our IILM points. The p meson contribution is needed to explain the nearly flat 
behavior of the effective mass for r > 0.6 fm. 



TABLE I: 7r, p, a± masses (in GeV units) calculated in the IILM for different quark masses. 



m q 




M p 


M ai 




r B1 


0.02 


0.30 ±0.04 


1.0 ±0.1 


1.6 ±0.1 


~ 0.01 


< 0.02 


0.03 


0.36 ± 0.04 


0.9 ±0.1 


1.6 ±0.1 


i$ 0.01 


< 0.03 


0.05 


0.46 ± 0.04 


1.0 ±0.1 


1.7 ±0.1 


~ 0.05 


< 0.01 


0.07 


0.53 ±0.04 


1.0 ±0.1 


1.7 ±0.1 


~ 0.05 


< 0.01 


0.09 


0.60 ± 0.04 


0.9 ±0.1 


1.8 ±0.2 


~ 0.05 


< 0.01 
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FIG. 3: Effective mass plots in the p meson channel evaluated in the IILM at different quark masses and compared to the 
phenomenological parametrization (solid line). The dashed straight line represents the expected asymptotic plateau if the p 
meson decays into two pion, while the solid straight line represents the expected asymptotic plateau if the p meson decays into 
two constituent quarks with masses estimated from the ILM calculation in [23| . 



Our results for the calculation of the a\ meson effective mass for several quark masses are shown in FigHJ where 
they are compared with the best fit obtained from the phenomenological representation of the spectral function ^ . 
In order to reduce the number of free fitting parameters, we have chosen to neglect the contribution of the continuum, 
and to restrict the fit of the IILM points to the region r > 0.6 fm. In addition, we have used the value for the pion 
mass and decay constant calculated in the IILM in [n| . The qualitative behavior predicted in the previous session is 
very well reproduced by our model. In particular, we observe that the singularity arising from the cancellation of the 
7r and a\ contribution to the two-point correlator is clearly developed. This result provides a clean evidence that both 
the pion and the aq meson exist in the instanton vacuum. We note that, in this channel, the simple parametrization 
© of the ALEPH Collaboration data [2l[ is quite poor in the low s region. This is presumably the source of the 
small discrepancy observed for some quark masses, at the largest Euclidean times. 

For comparison, in Fig[5]we show our effective mass plots for the pion, which were obtained in T9| and display a 
completely flat behavior at large Euclidean times. 

The complete list of p meson and a\ masses extracted from the fit of the effective mass plot are summarized in 
Table HI In general the masses obtained in the IILM are about 30% larger than the corresponding experimental values. 
This discrepancy suggests that about 1 /4 of the resonance mass is due to correlations which are not related to chiral 
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FIG. 4: Effective mass plots in the a± meson channel, evaluated in the IILM at different quark masses and compared to the 
phenomenological parametrization (solid line). The latter does not include the contribution from the perturbative continuum 
(see discussion in the text). 



symmetry breaking. This hypothesis is confirmed by the fact that the chiral asymmetry, 

x= m ai -m p 
m ai + m p 

which has been suggested as a parameter quantifying the contribution of chiral forces to hadron mass splittings [24| , 
is remarkably well reproduced in this model. At the lightest quark mass we find y^ 1LM ~ 0.2, which should be 
confronted with the experimental value x exp ' — 0.23. 

Let us now look in more detail on the mechanism of meson decay and we focus on the p meson in the IILM. The 
effective mass plots in the vector meson channel calculated at different quark masses are consistent with a small width 
spanning from 10 MeV to 50 MeV. Interestingly, the vector meson is always unstable in this model, even at quark 
masses for which there is no phase-space for decaying into pions, since 2M„ > M p . A natural explanation of this fact 
is that, in the absence of confinement, these hadrons decay into their quark-antiquark constituents. In the instanton 
vacuum, quarks acquire an effective mass as a consequence of the spontaneous breaking of chiral symmetry. Such a 
"constituent" quark mass at rest M q was calculated in several approaches and found to range from ~ 400 MeV to 
200 MeV, depending on the bare quark mass (see e.g. (2^] and references therein). In Fig. [31 we compare the effective 
mass plot for the p meson with the 2 M n line, and the 2M q line calculated in our model. We can see that, for the 
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FIG. 5: EfFective-mass plots for the pion evaluated in the IILM in [ll| 



two largest quark masses, the effective mass falls below the 2M W line, but always remains above the 2M q line. Hence, 
there is always phase space available to decay into constituent quarks. 



IV. CONCLUSIONS 



In this work we have studied the contribution of instanton-induced chiral forces in the a\ and p meson resonances. 

We have provided clean evidence that the a\ and p meson can exist even in the presence of instanton forces only 
(i.e. in the instanton vacuum). Their masses extracted from a fit of the effective mass plot have been found to be 
about 30% larger than the corresponding experimental values. This deviation should be confronted with the excellent 
agreement with the experimental masses previously found in the case of the pion and the nucleon (l9l | . These results 
suggest that chiral forces are weaker in these resonances, but still represent the leading source of interactions. 

On the other hand, the confining forces are very important to determine the width of these resonances and their 
decay properties. In fact, we have observed that, in this model, the vector meson is unstable even when quark masses 
are large and there is no phase space to decay into two pions. This can be explained assuming that — as a consequence 
of the absence of confinement — mesons can dissociate into their constituents. 

In the future, it would be interesting to study how the stability of the meson is restored, once confinement is 
introduced in the model. A way to do so would be to extend the ensemble of gauge configurations to include 
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regular-gauge instantons and to study the behavior of the effective mass as a function of the density of such pseudo- 
particles [251 ] . 
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